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Abstract 

, We analyze the exclusive rare B -^ K^*'i^£^, B -^ K^*'vl> and B -^ K*j decays in the 

Ph. 

(— I I Applequist-Cheng-Dobrescu model, which is an extension of the Standard Model in presence of 



universal extra dimensions. In particular, we consider the case of a single universal extra dimen- 
sion. We study how spectra, branching fractions and asymmetries depend on the compactification 
parameter 1/R, and whether the hadronic uncertainty due to the form factors obscures or not such 
a dependence. We find that, using present data, it is possible in many cases to put a sensible 
lower bound to 1/R, the most stringent one coming from B -^ K*^. We also study how improved 
experimental data can be used to establish stronger constraints to this model. 

PACS numbers: 12.60.-1, 13.25.Hw 



I. INTRODUCTION 

Although the Standard Model (SM) of electroweak interactions has successfully passed 
several experimental tests, it is commonly believed that a more fundamental theory should 
exist. Direct evidence of new Physics will be hopefully gained at high energy colliders such 
as the LHC In the meanwhile, signals of new interactions and particles can be obtained 
indirectly through the analysis of processes which are rare or even forbidden in the Standard 
Model. Among these, rare B decays induced hj b -^ s transition play a peculiar role since 
they are induced at loop level and hence they are suppressed in the Standard Model. 

Present data show that such a suppression indeed occurs. In the case of 6 — *> 57 modes, 
which in SM are induced by the electromagnetic penguin diagrams dominated by top and 
W exchange, the branching fractions have been measured both for inclusive and exclusive 
transitions; we collect them in Table HI 

TABLE I: Branching fractions of rare B decays induced by 6 ^ 57 transition. 



Mode Belle Collab. 


BaBar Collab. 




B -^ X,7 (3.55 ± 0.32 igj? ±[J:J^) x 10" 


Ml] (3.27 ±0.18 igi^ ±0.04) ^10- 


'[2] 


B°^K*°-f (4.01 ± 0.21 ± 0.17) X 10-5 


[3] (3.92 ± 0.20 ± 0.24) X 10-5 


m 


B- ^ K*--f (4.25 ± 0.31 ± 0.24) X 10-5 


[3] (3.87 ± 0.28 ± 0.26) X 10-5 


[4] 



Modes with two leptons in the final state, such as B ^ Xgi^i^, B -^ K^*H'^i^ and 
B -^ Xsi'i', B -^ i^*-*Vz/, are also suppressed, b -^ si'^i" transitions are described in SM 
by QCD, magnetic and electroweak semileptonic penguin operators, which give rise to an 
effective Hamiltonian composed of ten operators, as we shall see in more detail below. The 
resulting SM predictions depend on both the chirality structure of such operators, both on 
the value of their Wilson coefficients. The situation is simpler in the case oib ^ svv modes, 
described by Z penguin and box diagram dominated by top exchange: the corresponding 
effective Hamiltonian is composed by a single operator, therefore new Physics effects can 
just induce an operator with opposite chirality or modify the value of the Wilson coefficient, 
a scenario relatively simple to analyse. 



From the experimental point o 



view, the most recent measurements of the branching 



fractions are provided by Belle p, |7| and BaBar ja |8| Collaborations and are collected 
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TABLE II: Branching fractions of rare B decays induced by 5 ^ sl'^i and b — > sw transitions. 



Mode 


Belle Collab. 






BaBar Collab. 




B -^ Xsi+i- 


(4.11±0.83±[]|f) X 10-6 


m 


(5, 


,6 ±1.5 ±0.6 ±1.1) X 10-6 


[6J 


B -^ K*e+£- 


(16.5 ±|i ±0.9 ±0.4) X 10-7 


m 




(7.8 ±}:? ±1.2) X 10-7 


[8] 


B -^ Ki+i- 


(5.50 ±lfQ ±0.27 ± 0.02) x 10-' 


^ [7] 




(3.4 ±0.7 ±0.3) X 10-7 


m 


B -^ Xsuv 


B -^ K*vv 


B- -^ K-u9 


< 3.6 X 10-5 (90% CL) 


f9] 




< 5.2 X 10-5 (90% CL) 


[10] 



in Table |n] In addition, important information can be gained from the forward-backward 
lepton asymmetry m B ^ K*£^i^ , which is a powerful tool to distinguish between SM and 
several extensions of it. Belle Collaboration has recently provided the first measurement of 
such an asymmetry |ll| . 

Among the various models of Physics beyond the SM, those with extra dimensions are 
attracting interest for manifold reasons. For example, they provide a unified framework for 
gravity and other interactions, giving some hints on the hierarchy problem and a connection 
with string theory. Particularly interesting are the scenarios with universal extra dimensions, 
in which all the SM fields are allowed to propagate in the extra dimensions. Their feature 
is that compactification of the extra dimensions leads to the appearance of Kaluza-Klein 
(KK) partners of the SM fields in the four dimensional description of the higher dimensional 
theory, together with KK modes without corresponding SM partners. A simple scenario is 
represented by the Appelquist-Cheng-Dobrescu (ACD) model [12j with a single compactified 
extra dimension, which presents the appealing feature of introducing only one additional free 
parameter with respect to the Standard Model, i.e. 1/R, the inverse of the compactification 
radius. 

Analyses aimed at identifying the signatures of extra dimensions in processes already 
accessible at particle accelerators or within the reach of future facilities give different bounds 
to the sizes of extra dimensions, depending on the specific model considered. The bounds 
are more severe in the case of UED, and in the 5-d scenario constraints from Tevatron run 
I allow to put the bound 1/R > 300 GeV. In the following we analyze a broader range 



\/R > 200 GeV to be more general. 

n 

Rare B transitions can be used to constrain the ACD scenario [13[. In particular, Buras 
and collaborators have investigated the impact of universal extra dimensions on the B^^ — 
B^ ^ mixing mass differences, on the CKM unitarity trian gle and on inclusive b —>■ s decays for 
which they have computed the effective Hamiltonian J14l.ll5l|. The availability of precise data 
on exclusive b -^ s- induced decays, collected in Tables IHTTl induces us to extend the analysis 
to such processes in the framework of the ACD model. In this case, the uncertainty in the 
hadronic form factors must be taken into account, since it can overshadow the sensitivity to 
the compactification parameter. Actually, we show that this is not the case, at least for the 
smallest values of l/R: computing, for example, the branching fractions of 5 — > K'^*^i^i~ 
as well as the forward-backward lepton asymmetry in i? -^ K*i^i^ for a representative 
set of form factors we find that a bound can be put, and it can be improved following 
the improvements in the experimental data. Moreover, since in the limit of large l/R 
the Standard Model is recovered, we also investigate the agreement of current data with 
SM predictions [16]. We have also considered the decays B -^ K^*^ui', although for these 
modes no signal has been observed, so far, studying how observables like the various helicity 
amplitudes for B -^ K* transitions are modified in the ACD model. Finally, we have 
considered the branching fraction oi B -^ K*'y as a function of l/R, pointing out that it 
allows to establish the most stringent bound for the compactification parameter. 

The plan of the paper is as follows: in the next Section we briefly describe the ACD 
model. We discuss the modes B -^ K^*H^i~ , B -^ K^*^ui' and B -^ K*'y in the subsequent 
Sections; finally, we present our conclusions and the perspectives for further analyses. 

II. MODELS WITH EXTRA DIMENSIONS: THE ACD MODEL WITH A SIN- 
GLE UED 

If other dimensions exist in our universe apart the usual 3 spatial + 1 temporal ones, 
and if such extra dimensions are compactified, fields living in all dimensions would mani- 
fest themselves in the 3+1 space by the appearence of Kaluza-Klein excitations (from the 
original Kaluza and Klein studies aimed at unifying electromagnetism and gravity by the in- 
troduction of one extra dimension [13]). For example, in the case of a single extra dimension 
with coordinate x^ = y compactified on a circle of radius R (the compactification radius), a 



field F{x,y) (with x denoting the whole of the usual 3+1 coordinates) would be a periodic 
function of y, and hence it could be expressed as 

n=+oo 

F{x,y)= Y. i^n(x)e^"-^/^ . (2.1) 

n=— oo 

If, for example, F is a boson field obeying the equation of motion OmQ^^ F{x,y) = 
(M=0,l,2,3,5), the KK modes Fn would obey the equation 

(d^d, + ^] F^{x) =0 /i = 0, 1, 2, 3 , (2.2) 

and therefore, apart the zero mode, they would behave in four dimensions as massive particles 
with ml = [n/RY. 

An important question is whether ordinary fields propagate or not in all the dimensions. 
One possibility is that only gravity propagates in the whole ordinary + extra dimensional 
Universe, the "bulk". Opposed to this are models with universal extra, dimensions (UED), 
in which all the fields propagate in all available dimensions. 

In this paper we focus on the model developed by Appelquist, Cheng and Dobrescu 
(ACD) [l2] that belongs to the UED scenarios. It consists in the minimal extension of the 
SM in 4 + (5 dimensions, and we consider the simplest case 6 = 1. This extra dimension is 
compactified to the orbifold S^/Z2, with the coordinate x^ = y running from to 2ttR. The 
points y = 0,y = ttR are fixed points of the orbifold; the boundary conditions at these points 
determine the Kaluza Klein mode expansion of the fields. Under the parity transformation 
P5 '■ y ~^ "V fields having a correspondent in the 4-d SM should be even, so that their zero 
mode in the KK mode expansion can be interpreted as the ordinary SM field. On the other 
hand, fields having no correspondent in the SM should be odd, and therefore they do not 
have zero modes in the KK expansion. For example, in this scenario a vector field has a 
fifth component which is odd under P5. 

Important features of the ACD model are: i) there is a single additional free parameter 
with respect to the SM, the compactification radius R; ii) the conservation of KK parity, 
which has the consequence that there is no tree-level contribution of KK modes in low energy 
processes (at scales /i ^ 1/-R) and no production of single KK excitation in ordinary particle 
interactions. 

A detailed description of the extension of SM in five dimensions is provided in [ij]; here 
we recall the main features of such a construction. 



Gauge group 

The gauge bosons associated to SU(2)l x U{1)y gauge group are W^ (a = 1,2,3, 
M = 0, 1, 2, 3, 5) and Bm, and the gauge couphngs are: ^2 = fi'2V27ri? and g' = g' \J'1t\R 
(we denote with the caret the quantities referring to the five dimensional description). 
The charged bosons are W^ = -^(Wlj =F W^j). Moreover, as in SM, W^f and Bm mix 
giving rise to the fields Zm and Am- The mixing angle is defined through the ordinary 
relations: 

cw = cos^H^ = . sw = sin 611^ = . . (2.3) 

Due to the relations between the five and four dimensional gauge constants, the Wein- 
berg angle is the same as in SM. On the other hand, the gauge bosons associated to 
SU{?>)c are the gluons G'lj{x, y) (a = 1, .., 8). 

Higgs sector 

The Higgs doublet is written as: 

'^^ (2.4) 

75(^-^^3); 

where x^ = 75 (-^^ "F -^^)' ^'^o'^g these fields only tp has a zero mode, and we assign 
to such a mode a vacuum expectation value v, so that ip -^ v + H. H can be identified 
with the SM Higgs field, while the vacuum expectation value in 5 dimensions is related 



to the corresponding quantity in 4 dimensions through the relation: v = v/V^nR. 

Mixing between Higgs fields and gauge bosons 

The charged W^^^ and xfn) fields mix, as well as the neutral ^5(„) and X3(n)- The 
resulting fields are G9 -, , Gr^^-j which are Goldstone modes giving mass to the W^^Jt and 
Z,'^-), and a9 -), a^^, new physical scalars. 

Yukawa terms 

In order to have chiral fermions as in the SM, the left and the right-handed components 
of a given spinor cannot be simultaneously even under P5. The Yukawa coupling of the 
Higgs field to the fermions provides the fermion mass terms, and the diagonalization 
of such terms leads to the introduction of the CKM matrix, with the same steps as 



in the SM. In this respect, the ACD model belongs to the class of minimal flavour 
violation models, since there are no new operators beyond those present in the SM and 
no new phases beyond the CKM phase. As a consequence, Buras and collaborators 
have shown that the unitarity triangle is the same as in the SM [14]. In order to obtain 
4-d mass eigenstates for the higher KK levels, a further mixing is introduced among 
the left-handed doublet and the right-handed singlet for each flavour /. The mixing 

TTl f / Tl 

angle is such that tan(2a /•(„)) = — ■/- {n > 1) giving the masses mftn) = \ m'j + — r, 

n/R V R 

so that it is negligible for all flavours except the top. 

As a result of the construction, the four-dimensional Lagrangian, obtained integrating 
over the 5th dimension y: 

C,{x)= r"" C,{x,y) (2.5) 



/o 
describes i) zero modes (corresponding to the Standard Model fields), ii) their massive KK 

excitations, iii) KK excitations without zero mod es ( they do not correspond to any field in 

SM). The related Feynman rules can be found in |14] . 



III. DECAYS B -^ K(*k+i- 

In the Standard Model the effective AB = — 1, AS* = 1 Hamiltonian governing the rare 
transition b -^ si'^i~ can be written in terms of a set of local operators: 

G ^^ 

Hw = 4 -^VttV,: Y: a(/^)0.(/i) (3.1) 

where Gp is the Fermi constant and Vij are elements of the Cabibbo-Kobayashi-Maskawa 

V V* 
mixing matrix; we neglect terms proportional to VubV*^ since the ratio 



is of the 



VtbVtl 
order 10~^. The operators Oj, written in terms of quark and gluon fields, read as follows: 

Oi = {sLaY^La){cLI3lf,CLf5) 

03 = {sLa'^r^bLa)[{uLf3'^r^lULf3) + ■■■ + (bL|3l^lbL|3)] 

04 = {sLa1^bLi3)[{uL(^-ff,ULa) + ■■■ + (bLl3j^bLa)] 

05 = isLal''ha)[{uRf3'y^,URp) + ... + {bnf3lt,bnf3)] 
Oq = {sLal^bLl3)[{uRp-ff,URa) + ... + (bR^-f^bRa)] 



^^7 = T?r^^b{sLaCr^''hRa)F^^ 



2 /a/3"^'^ 






IQlT 

^10 = TF^(^"i«7'^feLa) ^^7^75^ (3.2) 

1 ih ^^ i 

where a, f3 are colour indices, bjiL = b, and o"^*^ = -[7^,7*^]; e and (/^ are the 

electromagnetic and the strong coupling constant, respectively, while F^j/ and CJ^^ in O7 and 
Og denote the electromagnetic and the gluonic field strength tensor. Oi and O2 are current- 
current operators, O3, ..., Oq are QCD penguin operators, O7 (inducing the radiative b ^ S'j 
decay) and Og are magnetic penguin operators, O9 and Oio are semileptonic electroweak 
penguin operators. 

The Wilson coefficients appearing in ()3.ip have been computed at NNLO in the Standard 
Model [l9| . At NLO the coefficients have been computed also for the ACD model including 

nn 

the effects of KK modes J14l . I15| : we use these results in our study. No operators other 
than those collected in eq. ()3.2|l are found in ACD, therefore the whole contribution of the 
plethora of states only produces a modification of the Wilson coefficients that now depend 
on the additional ACD parameter, the compactification radius. For large values of 1/R the 
Standard Model phenomenology should be recovered, since the new states, being more and 
more massive, decouple from the low-energy theory. Our aim is to establish a lower bound 
on 1/R from the various B -^ K^*H^£^ observables. 

In the following we do not consider the contribution to -B ^ K*i^i~ with the lepton 
pair coming from cc resonances, which is mainly due to the operators Oi and O2 in ()3.2p . 
It can be experimentally removed applying appropriate kinematical cuts around the reso- 
nances. QCD penguins 03,...,0q can also be neglected since their Wilson coefficients are 
very small com par ed to the others. Therefore, we only need the coefficients C7, Cg and Ciq: 
as discussed in (ij, |l^ , the impact of the KK modes consists in the enhancement of Cio and 
the suppression of C7. 

The Wilson coefficients in ACD are modified because particles not present in SM can 

contribute as intermediate states in penguin and box diagrams. As a consequence, the 

2 

Wilson coefficients can be expressed in terms of functions F{xt,l/R), Xt = 1 , which 



generalize the corresponding SM functions -Fo(xt) according to: 

oo 

F{xu l/R) = Fo{xt) + Y. Fn{xu Xn) (3.3) 



n=l 



mt , n 



where Xn = 2 ^^^ "^^ ~ 15- The relevant functions are the following: C{xt, 1/-R) from Z^ 



i-W 



penguins; D{xt-, 1/-R) from 7 penguins; -E'(xj, 1/-R) from gluon penguins; D'{xt, l/R) from 7 
magnetic penguins; E'{xt, l/R) from chromomagnetic penguins. The functions can be found 



m 



m 



15|; here we collect the formulae needed in our analysis. 



C7 

In place of C7, one defines an effective coefficient C7 which is renormalization 

scheme independent |l8| : 

C?''"{f.,) = V^^Cr\f^w) + I (^- - r/i) C^'\f^w) + Ci'\^w) E M- (3.4) 

as(/iiy) 

where 77 = — - — — , and 
as(/i6) 

Ci'\fiw) = 1 C^'\fiw) = -\d\xu l/R), C^ifiw) = -lE'ixt, l/R); (3.5) 

the superscript (0) stays for leading log approximation. Furthermore: 





Ol = 


14 
23 


16 

"^ = 23 


6 

a3=23 "^ = 


12 
23 


as 


= 0.4086 


ae = 


-0.4230 


aj = -0.8994 


as = 0.1456 


h 


= 2.2996 


h2 = 


-1.0880 


h = -ij hi-- 


1 
14 


h = 


-0.6494 


he = 


-0.0380 


hj = -0.0185 


/ig = -0.00. 



(3.6) 



The functions D' and E' are given by eq. (j3.3j) with 

^°^'^*^- 12(1 -a;,)^ + 2(l-x,)4 ^'''^* ^^'^^ 



D'niXuXn) 



Xt{-37 + AAxt + llxl + 6x2 (10 _ g^^ ^ ^^2-^ _ 33,^(21 - 54xf + 17x?)) 



36(xt - 1)3 



+ 



Xn\^ iXn + oXj^) Xr, 



6 1 + Xn 

-2 + Xn + 3Xt) {Xt + SX^ +xl{3 + Xt) - Xn{l + (-10 + Xt)Xt)) Xn + Xt 



In- 



^n\Xtl Xf; 



6(Xi - 1)4 

xj(-17 - Sxt + x1 - 3xn(21 - Qxt + x^) - Gx^ (10 - 9xt + Sx?)) 
~ 12(xt - 1)3 

- -x„(l + x„)(-l + 3x„)ln- 

2 l+Xn 

[l + x„)(xt + 3x^2 + x^(3 + Xt) - x„(l + (-10 + Xi)xi)) ,^ x^ + Xf 

l + Xn 



l + Xn 

(3.9) 



+ 



M 



2(xt - 1)4 



In- 



.(3.10) 



Following [14| one gets the expressions for the sum over n: 



2^ Dn{Xt,Xn) 



n=l 



nMwR 



Xt(-37 + Xi(44 + 17xf)) 
72(xi - 1)3 

1 (2|/V2 + 7^3/2 ^ 3^5/2) 



dV' 



+ 



6 

-2 + 3x^x^(1 + 3x0 
Q{xt - 1)4 



coth(7rMH/i?v^) 



J(i?,-l/2) 



[Xi(l + 3xi) - (-2 + 3xi)(l + (-10 + xt)xt)] J{R, 1/2) 



6(Xi - 1)4 

[(-2 + 3xi)(3 + Xt) - (1 + (-10 + xt)xt)] J{R, 3/2) 



6(xi - 1) 

(3 + Xt) 



-J(i?, 5/2) 



E K(^*, 



Xr, 



n=l 



+ 



ttMwR 



6(Xi-l)4' 

Xt(-17+(-8 + Xt)xt) 
24(xi - 1)3 

,1/2 , ^^3/2 



(3.11) 







rfy (y^/2 + 2|/3/^ - 3y5/2) coth(7rM^i?v^) 



Xf(l + 3xt) 
(x,-l)4 ^(^'-V2) 

J [xt(l + 3xi) - (1 + (-10 + xt)xt)] J{R, 1/2) 



+ 



(xt - 1) 
1 

{Xt - 1)4 

(3 + X,) 

(X* - 1)4 



[(3 + Xt) - (1 + (-10 + xt)xt)] J{R, 3/2) 
J(i?, 5/2) 



(3.12) 



where 



J{R,a)= / dyy'^ coth{Tr MwRy/y) 



X 



l+o 



coth^nTntRy/y) 



(3.13) 



In the ACD model and in the NDR scheme one has 

CM = P^'''' + ^^"^^^y^^ - 4Z(x„ l/R) + P^E(x„ l/R) (3.14) 

sm t/vK 

where Pq^^ = 2.60 ± 0.25 [20| and the last term is numerically negligible. Besides 



with 



and 



Yixt,l/R) = Yoixt) + J2Cnixt, 



Xr. 



n=l 

oo 



Z{xt, l/R) = Zo{Xt) + Y^ Cn{Xt,Xn) 

n=l 



Yoixt) = ^ 



Xi-4 



3xt 



Xt 



Zo{xt) 



1 {Xt 



\\iXt 



18x^ - 163x3 + 259x2 - lOSx* 



+ 



144(xi - 1)3 
32x^ - 38x3 - 15x2 + 18xi ^ 
72(xi - 1)4 9 



Inxt 



Xt 



Ka^t - 1)' 



X, 



^Xf + 7 + (3 + 3xf + 7x„ - XfX„) In 



Xt + Xr. 
l+Xn 



(3.15) 



(3.16) 
(3.17) 



7 , Cn[Xt, 



Xr. 



n=l 



Xt(7 — Xt) nMwRxt 



I6{xt - 1) 16(xt - 1) 



[3(l + Xt)J(i?,-l/2) + (xi-7)J(i?,l/2)] . (3.18) 



C 



10 



Cio is fj, independent and is given by 



C 



Yixt, l/R) 



10 



sin^^ 



w 



(3.19) 



We fix the renormalization scale to /i = /if, ~ 5 GeV. 

With these coefficients and the operators in ()3.2|) the inclusive b — > si^i~ transitions have 
been studied in [ij, |l5|. The exclusive B —^ K^*'^i'^£~ modes, on the other hand, involve 



10 



the matrix elements of the operators in ()3.2|) between the B and K or K* mesons, for which 
we use the standard parametrization in terms of form factors: 



< Kip')\rs^,b\Bip) >= {p + p'),F,{q^) + ^^ ,^^ g, (Fo(g^) - F,{q^ 



g = p-p', Fi(0)=Fo(0))and 

< K{p')\-sia,,q^h\B{p) >= \ip + p'),q' - (M^ - M'^)q, 



FAq' 



Mb + Mk' 



(3.20) 



(3.21) 



<K*{p',e)\sj,{l-j,)b\B{p)> = e^,„^e*>V^ 



— I 



Mb + Mk* 
el{MB + MK^)A,{q')-{e*-q)ip + p'] 



A2{q' 



2M 



K 



and 



[e* ■ q)^{A,{q') - Ao{q'))q, 



<K*{p',t)\m^,q^^-^^-^b\B{p)> = le^^^pe^'^py'^ 2T,{q') + 



^{Mb + Mk*) 
(3.22) 



+ 



e;(M|-M|0-(e*-g)(p + p') 



n(Q^) 



+ (e* • q) 



A^ can be written as a combination of Ai and A2: 



% 



Ml - Ml, 



-{p + p') 



n{q') . 



(3.23) 



^3(g ) = ^777 Mq ) 777 Mq ) 



2Mk* 



2Mk' 



(3.24) 



with the condition ^3(0) = Ao(0). The identity a^ul^ = — ^e^i/a/jcr"^ (^0123 = +1) iniphes 
that Ti(0) = r2(0). 

The form factors are non-perturbative quantities. We use for them two sets of results: 
the first one, denoted as set A, obtained by three-point QCD sum rules based on the short- 
distance expansion '21|; the second one, denoted as set B, obtained by QCD sum rules based 
on the light-cone expansion J22]. In both cases, the mass of the h quark is finite. The q^ 
dependence is fitted in the region where the methods can be relibly applied, actually the 
low q^ region, and then it is extrapolated to the full physical range. The form factors in set 



11 



A are fitted with by two different functional dependences: either polar or linear. Fi, V and 
Ti display a polar behavior, while Ai, T2, and A2, T3 depend linearly on g^, with decreasing 
(increasing) behaviour. Only the form factor Ft is a double pole. The values of parameters 
with their estimated errors can be found in J2ll]; the uncertainties have been included in the 
analysis we present below. In set B of form factors, Ai and T2 are fitted as simple poles, 
V and Ti as a sum of two polar functions and Fi, Ft, A2, T3 are the sum of a pole and a 
double pole. The values of the parameters and of their estimated uncertainties can be found 



m 



22| : also for this set we include in the numerical analysis the errors of the parameters. 



The main differences between the two sets of form factors essentially concern Ai and To. 

n 

In the numerical analyses we use the values reported by the PDG j23| for masses and 
CKM matrix elements. We also use mh = 4.8 GeV, rrit = 172.7 GeV, coinciding with 



the central value recently reported by the Tevatron electroweak working group 
tbo = 1.527 ± 0.008 ps ^. 



2J, and 



B -^ K£+t 



For each set of form factors the differential decay rate in the invariant mass squared of 
the lepton pair 



dr 

dq 



.(B 



") 



MlGla 




(3.25) 



(g^ = M|^^_) displays a dependence on the compactification parameter l/R, as depicted 
in Fig. n where we have considered the values l/R = 200, 500 GeV and the case of the 
Standard Model (large values of l/R)- The maximum effect is in the low g^ region, where the 
spectrum has the maximum. However, such an effect is obscured by the hadronic uncertainty, 
which is evaluated comparing the two set of form factors and taking into account their errors. 
Theferore, the differential decay rate does not seem the most suitable observable for studying 
the effect of extra-dimension at the current level of hadronic uncertainties. The situation is 
different for the width. In Fig. |21we plot, for the two sets of form factors, the branching 



12 



fraction as a function of 1/R and compare it with the experimental data provided by BaBar 
and Belle . A constraint cannot be put on 1/R if one adopts set A, while set B allows to 
exclude 1/R < 200 GeV. It is interesting to observe that, within the Standard Model, set 
A prefers the lowest experimental range, corresponding to the BaBar result, while set B is 
in better agreement with Belle data. Improved measurements should resolve the present 
discrepancy between the two experiments. At the same time, they should reather easily 
allow to increase the lower bound for the compactification parameter. 




q^ (GeV^) 




q^ (GeV^) 



FIG. 1: Differential branching fraction dBR{B -^ Ki i^)/dq obtained using set A (left) and B 
(right) of form factors. The dark (blue) region is obtained in SM; the intermediate (red) one for 
1/R = 500 GeV, the light (yellow) one for 1/R = 200 GeV. The contribution of cc resonances is 
not displayed. 



B. B 



K' 



A great deal of information can be obtained from the mode B -^ K*i~^i~ by investigating, 
together with the lepton invariant mass distribution, also the forward-backward asymmetry 
Afb in the dilepton angular distribution, which may reveal effects beyond the Standard 
Model that could not be seen in the analysis of the decay rate. In particular, in SM, due to 
the opposite sign of the Wilson coefficients Cj and Cg, Ajh has a zero the position of which 
is almost independent of the model for the form factors 26[. 

Let us define 6£ as the angle between the i~^ direction and the B direction in the rest frame 
of the lepton pair (we consider massless leptons). The decay amplitude can be written as 
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FIG. 2: BR{B -^ K£^i~) versus — using set A (left) and B (right) of form factors. The two 

R 
horizontal regions correspond to the experimental data provided by BaBar (lower band) and Belle 

(upper band), see Table ITll 



sum of non interfering helicity amplitudes, the double differential decay rate reads as follows: 



dq^dcosOi 



13-7r5 



21% 



Ml 



{sin^9iAL+ 



+ q' 



;i + cose,f{Al + A^) + (1 - coseifiA^i + a^)] } (3.26) 



where A^ corresponds to a longitudinally polarized K* \ while ^+(_) represent the con- 
tribution from left (right) leptons and from K* with transverse polarization: e± = 

1 



A, 



Mh 



and 



B^{Ml - Ml, - q^) + B^X + D^{Ml - Ml, - q") + D^X \ (3.27) 

Ai = \X^'^{A-C)T{Bi-Di)\^ (3.28) 

Al = \X^'\A + C)t {Bi + /^i)P , (3.29) 

where A = A(M^, M|;,, g^). The terms A,C,Bi,Di contain Wilson coefficients and form 
factors: 



^=7'""^'(«') + ^»M, + Afe. 



c = c 



10 



V{q' 



Mb + Mk* 



Bi 



C-j 



Am,T2{q^){Ml - Ml,) + C^A^{q^){MB + Mk*] 



(3.30) 
(3.31) 
(3.32) 



14 



Bo = - 



Cr 



^rrib T2(g2) + g2 



n 



(Ml 



M? 



Co 



A2{q' 



K* 



Mb + Mk* 



Dr = CioA^{q^){MB + MK*] 



D, 



-C 



A2{q' 



10 



Mb + Mk' 



(3.33) 
(3.34) 
(3.35) 



The forward-backward asymmetry, defined as 



A'^iq' 



d^T 



dq^dcosO, 



-dcosOp — 



rf^r 



dq^dcosO, 



■dcosOp 



d'T 



-dcosOp + 



d'T 



(3.36) 



reads 



A'^{q') 



dq^ dcosOp J-i dq'^dcosO, 

3 2q\Al + A^ - A^ - A^) 



-dcosOp 



(3.37) 



AAl + 2q^{AL + Af + A^ + A^) ' 
We can now discuss the predictions of the ACD model for the branching ratio as well as 
for the lepton forward-backward asymmetry. The differential branching ratio is shown in 
Fig. 121 As in the case oi B ^ Ki^i^ . the spectrum is enhanced for lower values of l/R, 
however, due to the hadronic uncertainty, it is not possible to clearly disentangle the extra 
dimension effect. As for the total rate, depicted in in Fig. |3]for the two sets of form factors, 
set A does not allow to establish a lower bound on l/R, while for set B one gets again 
l/R > 200 GeV. The present discrepancy between BaBar and Belle measurements does not 
permit stronger statements, more precise data from both the experiments are expected. 
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FIG. 3: Differential branching fraction dBR{B -^ K*ti~)/d(^ using set A (left) and B (right) 
of form factors. The dark (blue) region is obtained in the SM, the intermediate (red) one for 
XjR = 500 GeV, the light (yellow) one for l/R = 200 GeV. 
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FIG. 4: BR{B -^ K*£^i ) versus — using set A (left) and B (right). The two horizontal regions 

-ft 
correspond to BaBar (lower band) and Belle (upper band) results, see Table HH 



More information conies from the forward-backward asymmetry. We show in Fig. the 
predictions for the SM, 1/R = 250 GeV and 1/R = 200 GeV. The zero of Afb is sensitive 
to the compactification parameter, so that its experimental determination would constrain 
1/R. 
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FIG. 5: Forward- backward lepton asymmetry in S — > K*i^£ versus -— using set A (left) and 

R 
B (right). The light (yellow) bands correspond to the SM results, the intermediate (red) band to 

1/R = 250 GeV, the dark (blue) one to 1/R = 200 GeV. 



We can elaborate on this point, since it is easy to see, using ()3.27|) - ()3.35|) in ()3.37|) . that 
the position of the zero of the asymmetry, sq, is determined by the equation: 



Re{C,) + ^C, \{Ms + MxO^ + {Ms - M^*)^ 
s [ \/(s) ^il-s) 



= . 



(3.38) 



16 



It is noticeable that, in the large energy limit of the final state light vector meson, a model 
independent prediction for the position of the zero of the asymmetry can be obtained. As 
a matter of fact, in this limit the form factors oi B ^ P,V obey spin symmetry relations 
, broken by hard gluon corrections to the weak vertex and hard spectator interactions. 



In the heavy-quark limit one can write 28| (see also 29|, |30|,|3l|): 



{P\^P Ti h\B) = Ci{E, fii) ^pifii, E) + T,(E, M, u;, /i„) ® 0^(u;, /in) ® 0p(m, /in) + . . . , (3.39) 

where F, is a generic Dirac structure and the dots stand for sub-leading terms in A/rrib. 
In the case of a vector meson V two functions ^_l.|| (depending on the Dirac structure Fj) 
appear in place of C,p- The matrix elements in ()3.39p get therefore two contributions. The 
first one contains the short-distance functions Ci, arising from integrating out hard modes: 
fij < nib, and a "soft" form factor ^p which does not depend on the Dirac structure of the 
decay current. The second term in ()3.39j) f actor izes into a hard-scattering kernel T, and 
the light-cone distribution amplitudes (ps and (pp. A still controversial question is to what 
extent the first contribution is numerically suppressed by Sudakov effects 3^, although it 
has been put forward that the second term in (13.391) should be subleading with respect to 

nfin 

the first one j30|, |33|, |3J| . Neglecting 0{as) effects the approximate symmetry relations 
mentioned above between the vector and tensor form factors for B ^ V transitions read 

mm 

T,{E) _ 1 Mb T,{E) _ {Mb + Mk^) 

V{E) 2Mb + Mk' Ai{E) 2Mb ' 

The use of eq. ()3.40p in ()3.38|) produces a form-factor independent result for the position of 
the zero of the asymmetry: using our numerical input parameters, one would obtain in the 
Standard Model Sq^^^ ~ 3.61 GeV^. On the other hand, taking into account corrections to 
the relations in ()3.4()|1 . the position of the zero of Aft moves to sq ~ 4.2 ± 0.6 GeV^ J36l |. 

The dependence of Sq on the compactification parameter is depicted in Fig. IHl for the 
sets A and B of form factors. Two considerations are in order. The first one is that the 
value of So in the Standard Model is only marginally consistent with the result obtained in 
J2a|, suggesting that further corrections could shift sq to smaller values. The second one 
concerns the sensible dependence of sq on the compactification parameter: in particular, the 
zero is pushed to low values by decreasing 1/R. Such a sensitivity in the exclusive channel 
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is analogous to the one observed in the inclusive B — > Xsi^i~ mode, and indicates that sq 
is particularly suited to constrain 1/R. 

At present, the analysis of the forward-backward lepton asymmetry performed by Belle 
Collaboration indicates that the relative sign of the Wilson coeffcients Cg and Cj is negative, 
confirming that Afb should have a zero jll|] . The accurate measurement of Sq in the exclusive 
B -^ K*i^£~ channel is therefore within the reach of current experiments. 
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FIG. 6: Zero of the forward-backward lepton asymmetry versus — obtained using the set A (left) 

R 
anf B (right) of form factors. In the plots the horizontal region represents the the value of sq 

derived in [Sal. 



IV. THE DECAYS B -^ K^*) 



uv 



As mentioned in the Introduction, among the various flavour changing neutral current- 
induced 6— quark decays the transitions induced by 6 — > suu play a peculiar role, both from 
a theoretical and an experimental point of view. Within the Standard Model these processes 
are governed by the effective Hamiltonian 



n 



a 



eff 



V2 zvrsm (t^iy) 



(4.1) 



obtained from Z^ penguin and box diagrams where the dominant contribution corresponds 
to a top quark intermediate state. In ()4.H) 6w the Weinberg angle. Ol represents the left-left 
four-fermiori operator Ol = 67'' (1 — 75)3 i^'y^{l — 75) z^. The function X has been computed 
in 37 1 and J38|, |39|; we put to unity the QCD factor 7]x 



38, 



39 



40[. Possible New Physics 
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(NP) effects can modify the SM value of the coefficient cl, or introduce one new right-right 
operator: 

Heff = ClOl + Cr Or (4.2) 

{Or = 67^(1 + 75)5 1^7^(1 + 75)^^), with Cr only receiving contribution from phenomena 
beyond SM |4l|. 

Another reason of interest for b -^ svv is the absence of long-distance contributions 
related to four-quark operators in the effective hamiltonian. In this respect, the transition 
to neutrinos represents a clean process even in comparison with the b —>■ sj decay, where 



long-distance contributions are expected to be present, although small 



42|. 



Within the Standard Model, form factors are needed to predict branching ratios and 
decay distributions for the exclusive modes B -^ K^*'^!? (see, e.g. (4lL l43l|). The inclusion 
of effects stemming from one universal extra dimension is straightforward and requires the 
generalization of the function X{xt) |l4| : 



X(Xt, l/R) = Xo(Xt) + Y. Cn{Xt, Xn) (4.3) 

n=l 

where: 



X,{x^) = I 
and Cn{xt,Xn) defined in eq. ()3.17p . 



Xt + 2 3xt - 6 , 
1- 7 TTT In Xt 

xt-i [xt - ly 



(4.4) 



A. B ^ Kvv 

It is interesting to consider the missing energy distribution in the decay B -^ Kvv. We 
define Emiss the energy of the neutrino pair in the B rest frame and consider the dimension- 
less variable x = Emiss/Ms, which varies in the range 

~^ <x <l-^/r (4.5) 



2 
with r = M'^/Mq. The differential decay rate is then given by 



dV{B^Kvv) \cL + CR\'\Fi{q')\ 



3 '^ ' "^'1^^^ ^' VA'(g2,M|,M|-), (4.6) 

where q^ = M^(2x — 1) + M|- and the sum over the three neutrino species is understood. 

In Fig. [7| we plot the SM missing energy distribution, together with the distributions 
obtained in ACD for different values of l/R. This distribution is sensitive to l/R, and the 
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FIG. 7: Missing energy distribution in the decay B -^ Kuu for set A (left) and B (right) of form 
factors. The sum over the three neutrino species is understood. The dark (blue) region is obtained 
within SM, the intermediate (red) one for 1/R = 500 GeV and the light (yellow) one for 1/ R = 200 
GeV. 
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FIG. 8: BR{B -^ Kuv) versus 1/R for set A (left) and B (right) of form factors. 



largest effect is in the low-x region, however the hadronic uncertainty is too large for envis- 
aging the possibility of constraining the compactification parameter. As for the branching 
fraction, its dependence on 1/R is shown in Fig. |H1 Only an experimental upper bound exist 
in this channel, which is presently too large for any consideration: however, considering Fig. 
IHlone sees that the 1/R dependence is too mild for distinguishing values above 1/R > 200 
GeV. 
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B. B ^ K*vv 



For this mode one can separately consider the missing energy distributions for longitudi- 
nally and transversely polarized K*: 

n2 



dx 



\CL? \P'\ 

247r3 M|-. 



{Mb + Mk*){MbE' - Ml.)A^ 



2M2 ,^„2 



Mb + Mk- 



wrMq' 



and 



dV± 
dx 



•'1^2 






2Mb|p'| 



-\/(g2)^(MB + M;,.)A(g' 



(4.7) 



(4.8) 



Mb + Mx* 
where p ' and £" are the K* three-momentum and energy in the B meson rest frame and 

the sum over the three neutrino species is understood. The missing energy distributions 

for polarized K* are shown in Fig. El while the unpolarized distribution and the branching 

fraction are plotted in Fig. [TUland[TTl respectively. 
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FIG. 9: Missing energy distribution in S — > K*vi> for a longitudinally polarized K* (lower curves) 
and a transversally polarized K* (upper curves) for set A (left) and B (right) of form factors. The 
sum over the three neutrino species is understood. The dark (blue) region corresponds to SM, the 
intermediate (red) one to 1/ R = 500 GeV and the light (yellow) one to \/R = 200 GeV. 



Both the polarized K* and total missing energy distributions depend on l/i?, but the 
hadronic uncertainty obscures such an effect as in 5 -^ Kvv. The dependence of the 
branching fraction is not strong for 1/ R > 200 GeV. 
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FIG. 10: Missing energy distribution for unpolarized K* , with the same notations as in Fig. |^ 
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FIG. 11: BR[B -^ K*^) versus 1/R using set A (left) and B (right) of form factors. 



V. THE DECAY B -^ K*j 

As a final case we consider the radiative channel B -^ K*'~f, which is the first one observed 
among such rare decay modes. The transition 6 -^ 57 is described by the operator O7 in 
the effective hamiltonian (j3.2p . and the B -^ K*'~f decay rate is given by: 

T{B -. K*j) = ^\V,,V:fml\ci'^'^^^^^m{0)rMl (1 - ^V . (5.1) 



Stt^ 



M| 



One can appreciate the consequences of the existence of a single universal extra dimension 
considering Fig. ^J where the branching fraction is plotted as a function of 1/R: the 
sensitivity to the compactification parameter is evident, and it allows to put a lower bound 
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of 1/ R > 250 GeV adopting set A, and a stronger bound of 1/R > 400 GeV for set B, which 
is the most stringent bound that can be currently put on this parameter from the set of B 
decay modes we have considered. 
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FIG. 12: BR{B -^ K*j) versus — using set A (left) and B (right) of form factors . The horizontal 

R 
band corresponds to the experimental result. 



VI. CONCLUSIONS AND PERSPECTIVES 

We have analysed the rare B -^ K^*^i~^i~, B -^ A'(*Vz/, B -^ K*-f decays in the ACD 
model with a single universal extra dimension, studying how the predictions for branching 
fractions, decay distributions and the lepton forward-backward asymmetry in B ^ K*i~^i~ 
are modified by the introduction of the fifth dimension. The possibility to constrain the only 
free parameter of the model, the inverse of the compactification radius 1/R, are slightly 
model dependent, in the sense that the constraints are different using different sets of 
form factors. Nevertheless, various distributions, together with the lepton forward-backward 
asymmetry in S — > K*i'^i~ are very promising in order to constrain 1/R. We found that 
the most stringent lower bound comes from B — > K*'-/. Improvements in the experimental 
data, expected in the near future, will allow to establish more stringent constraints for the 
compactification radius. 
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